Numerical study on discharge characteristics influenced by secondary electron emission in capacitive RF argon glow discharges by fluid modeling
1. IntroductionCapacitive RF glow discharges at low pressure have been widely adopted in many fields, such as microelectronics and material processing technology. Material processing technology with enhanced secondary electron emission (SEE) is adopted in various applications, such as electronic materials and nanomaterials. In order to improve the efficiencies of these applications, it is very important to understanding the generating mechanism and physical characteristics of the plasmas.
As is well known, the SEE is the emission of the electrons from the electrode surfaces, which occurs mainly when the accelerated ions collide with the electrodes. It leads to the SEE being one of the critical parameters in determining the magnitude of the glow discharge yield.[1–19] Its process under ion bombardment plays an important role in the plasma parameters such as electron density, electron temperature, electric field and potential, and ion and electron heating.[20–31] Many efforts have been made to explore the effect of SEE coefficient γ.[11–20] In the simulations and experiments of RF glow discharge studies, the effect of SEE coefficient γ on the sheath structure is very important, since the electron heating on the electrode surfaces is limited by the electric potential.[4] As the SEE increases, the sheath potential becomes smaller and its isolation effect of primary electron heating from the plasma to the electrode becomes weaker. The theoretical predictions concern the role of the electron heating at the sheath boundary in RF glow discharges, due to the oscillating sheath. The transition of the sheath structure can be controlled by increasing SEE.[4,5] The studies indicate that the sheath heating of the electrons could be the dominant heating mechanism in RF glow discharges, compared with the secondary electron heating.[23–25,31] Though these studies illustrate the discharge mechanisms with considering secondary electrons, the effect of SEE coefficient on the electron heating mechanism or on the discharge sustaining mechanism, which is a key question of RF glow discharge, has not been studied carefully.
In glow discharges, the mechanism of the electron heating largely depends on a wide range of discharge conditions such as electrode gap, gas pressure, driving frequency, applied voltage, etc.[18–23,29,30] Two mechanisms for the electron heating have been observed in the glow discharge: the electron pressure heating and the electron Ohmic heating. A few studies evaluated how effective the electron pressure heating is. The electron Ohmic heating is a continuous process, and the advantages of the electron Ohmic heating are as follows: i) the electrons can be heated homogeneously, and ii) the electrons with low thermal conductivity can also be heated. The SEE coefficient strongly affects ionization in glow discharge. In this paper, the effects of SEE coefficient on the electric field, electric potential, and electron temperature are studied. Also the SEE coefficient as a significant discharge parameter affects greatly on ion heating and electron heating, which consists of electron pressure heating and the electron ohmic heating.
In this work, a one-dimensional (1D) fluid model is used to study the effect of SEE coefficient on the RF glow discharge characteristics. The glow discharge is operated between two parallel plate electrodes. It is assumed that the electrode gap space is much smaller than the electrode diameter. Therefore, a 1D fluid model is reasonable.
The rest of this paper is organized as follows. In Section 2, the model is described. In Section 3, the results from the simulation are discussed. From our calculation results, we find that for all the parameters we choose in this work, after the 5000th cycle, the RF discharge reaches a steady state. Finally, some conclusions are drawn from the present study in Section 4.
2. ModelIn our model, a plasma of an RF glow discharge between two parallel-plate electrodes is described. When there exists a certain electric potential difference between the electrodes, the gas discharge between the electrodes will take place and the plasma will be generated. For a better study, we make several assumptions as follows.[31]
(i) The diameters of the electrodes are each much greater than the spacing between the two electrodes. Thus the discharge characteristics change only along the x axis.
(ii) The equations of continuity are written for the electron and ion density, the moment equations are obtained by the integrating Boltzmann equation over velocity space. The continuity equations are coupled to the Poisson equation relating the potential field to the particle density. The electron energy distribution is Maxwellian.
(iii) The diffusivity and mobility of ions and electrons can be described as continuum properties, i.e., continuity equations describe the creation, loss, and transport of charge carriers. Also both the ion flux and electron flux are described by the drift and diffusion approximation. By using a driving frequency of 13.56 MHz and a pressure of more than 100 mTorr (1 Torr = 1.33322×102 Pa), the momentum transfer collision frequency is more than the RF driving frequency. So the drift and diffusion approximation is valid.
(iv) Only the ground state ionization is considered, neither the metastable atoms nor any reactions with metastable atoms are taken into account.
(v) The temperatures of ions and neutral particles are constant and equal to gas temperature. Because of the inertia, ions cannot make a quick response to the change of electric field, i.e., ions can obtain little energy from the applied field.
The capacitively coupled RF argon glow discharge at low pressure is simulated. In the discharge, only the inelastic press for particles is considered:
.
Thus, the densities of the electrons and the ions satisfy the following equations:
where
and
are the fluxes of the electrons and the ions within the drift-diffusion approximation;
,
,
, and
are the diffusion coefficients and the mobilities of the electrons and the ions, respectively;
is the density of neutral atoms of gas and can be written as
(cm
−3), with
p being the pressure of neutral gas in units of Torr. The terms on the right-hand side in Eqs. (
1) and (
2) mean that in the source term, only the grounded ionization is considered. For the argon discharge, the ionization rate can be written as
where
is the electron temperature in units of eV. The electron energy balance is evaluated by
where
is the conductive heat flux of electron energy,
is the thermal conductive coefficient,
is the electron heating coming from the electron current coupled with the electric field as seen in Eq. (
9). The last term on the right-hand side in Eq. (
6) is the electron energy loss due to the ionization. The electron heating coming from the electron current has two parts: one is the electron heating coming from the electron pressure and the other is the electron heating coming from the ohmic heating, i.e.,
where
In this work,
is the electron pressure heating,
is the electron ohmic heating, and
is the electron heating. The electric field satisfies
where
e is the elementary charge and
ε0 is the vacuum permittivity of free space.
The initial conditions we used in this work are as follows:[28,31]
where
nε is the peak value of the initial charged particle densities,
ε is a small positive number, and
d is the gap between the parallel-plate electrodes.
The boundary conditions are as follows.[28,31]
at x = 0 (on the powered electrode):
,
,
, and
;
at x = d (on the grounded electrode):
,
,
, and
where γ and
are the SEE coefficient and electron recombination coefficient, respectively.
3. Results and discussionIn order to understand the plasma characteristics with SEE coefficients of 0.01, 0.1, 0.2, and 0.3 in CCRF glow discharge between parallel-plates electrodes, we numerically solve the model by a finite difference method with the boundary conditions and the initial conditions. The parameters used in the simulation are shown in Table 1.
Table 1.
Table 1.
![](https://cpb.iphy.ac.cn/html_resources/images/table-icon.gif) | Table 1.
Parameters used in the numerical simulation.
. |
Figure 1 shows the variations of cycle-averaged electric field, electron temperature, ionization rate, and electron density with SEE coefficient, at three positions: in the powered sheath at x = 0.02ċm, in the centre of discharge at x = 1.25 cm, and in the grounded sheath at x = 2.48 cm. Figure 1(a) shows that the electric field almost does not change with increasing SEE coefficients from 0.01 to 0.3. The electric field is virtually constant (zero) in the centre of discharge at x = 1.25 cm with all SEE coefficients. Figure 1(b) shows the cycle-averaged electron temperature changing with SEE coefficient increasing. It can be seen that the electron temperature in the centre of the discharge is almost constant for all SEE coefficients. While the electron temperature increases greatly in the two sheath regions compared with in the bulk plasma region as the SEE coefficient increases. This is because the secondary electrons gain high energy inside the sheath region, but these secondary electrons are unable to infiltrate to the bulk plasma region because the mean free path of electron is small under the present conditions.[31] Figure 1(c) shows the changes of the cycle-averaged ionization rate with SEE coefficient in a range from 0.01 to 0.3. The results show that the ionization rate has a significant change in discharge as the SEE coefficient increases; it has much higher change in the sheath region than in the bulk plasma region. The ionization rate is defined as the product of the electron density and the collision frequency between the electrons and neutral atoms. Although the electron densities in the two sheath regions are much smaller than in the bulk plasma region, but the electric fields in the two sheath regions are much larger than in the bulk plasma region. Therefore, as the SEE coefficient increases, the electrons in the two sheath regions can obtain more energies and the collision frequency between the electrons and the neutral atoms is very large, which leads to the larger ionization rates in the two sheath regions. Thus, the electrons before moving into the small ionization region, can generate more and more ionization rate in the sheath region.[16,18] Figure 1(d) shows that the higher electron density occurs in the centre of discharge as the SEE coefficient increases. The electric field slightly increases compared with the large increase of the electron density. The increasing electron density is consistent with the increase of electron temperature and the decrease of sheath thickness slightly as the SEE coefficient increases.[31]
Figure 2(a) gives the variations of cycle-averaged electron current density with x for the SEE coefficients of 0.01, 0.1, 0.2, and 0.3, respectively. It shows that the electron current density is positive in the powered sheath region and negative in the grounded sheath region, according to the electric field negative and positive values. The electron current density is very low in the center of the discharge because the electric field is very weak in the center of the discharge. The sign of the electric field changes when it goes from the powered electrode to the grounded electrode. The electron current density rapidly decreases at the location of low electric field in the sheath region. With the increasing of the SEE coefficient, the electron current density increases and has a peak in the powered sheath region. In the grounded sheath region, the electron current density also increases negatively and has a peak. This is due to the curvatures of the electric field lines in these regions. In the grounded sheath region, where the ion current density positively increases, so there are no ions emitted from the surface. It is clear that the electron current density plays a basic role in sustaining the RF discharge.
Figure 2(b) shows that the cycle-averaged ion current density negatively increases with the increasing of SEE coefficient in the powered sheath region. In the grounded sheath region, the ion current density increases positively as the SEE coefficient increases. It can be explained as follows. Owing to the high electric field and large electric field variations in the sheath region, the energy gained by the electrons from the electric field is not locally balanced by the energy loss through collisions, so that a significant electron flux of high-energy electrons can enter into the negative flow, the powered electrode acts as an electron gun. On the powered electrode, the fall energy gained from the electric field is much larger than the energy loss due to the collisions, so it is dominant in the negative glow.[23] The more the SEE coefficient increases, the more the ion flux created in the powered sheaths region is. The curve of ion current density has no peak-like electron current density because the ions are further accelerated as they pass through the sheath region.
The total current density is given by
. The change of the total current density is because the distribution of the electric field in the discharge is strongly affected in the presence of charge particles. But the electric field can also reach quite a large value in the range as indicated in Fig. 1(a) and does not depend on the total current density in the considered range as shown in Fig. 1(a). The cycle-averaged total current density is shown in Fig. 2(c). It should be noted that the total current density changes with SEE coefficient. The figure shows that the total current density is lower and negatively increases in the powered sheath region. In the grounded sheath region, the total current increases positively as SEE coefficient increases.
The secondary electrons flux on the surface of electrode is given by
. Figure 3 shows the cycle-averaged secondary electron fluxes as a function of the SEE coefficient. Firstly, at a position of x = 0.00 cm on the surface of the powered electrode, the electron current is positive and increases as the SEE coefficient increases, so more ohmic heating may take place there. At the same location, the ion current is negative and increases negatively as the SEE coefficient increases. Secondly, the secondary electron flux is negative on the surface of the grounded electrode, because the electron current increases negatively, while the ion current increases positively as the SEE coefficient increases. In the bulk plasma region the electron flux and the ion flux go through zero.[3,4,16]
The evolutions of the electron flux, secondary electron flux, and ion flux on the surface of the powered electrode in one cycle are shown in Fig. 4. From Fig. 4(a), the electron fluxes increase slowly from negative to zero in the time range from 0 t/T to 0.15 t/T, quickly from zero to positive in the time range from 0.15 t/T to 0.45 t/T; the electron fluxes decrease quickly from positive to zero in the time range from 0.45 t/T to 0.85 t/T, slowly from zero to negative in the time range from 0.85 t/T to 1 t/T. The evolution of the electron flux decreases with increasing SEE coefficient in the time ranges from 0 t/T to 0.15 t/T and from 0.815 t/T to 1.0 t/T. In the time range from 0.15 t/T to 0.815 t/T, the evolution of the electron flux increases largely as the SEE coefficient increases. It can be explained as follows: at the time of maximum sheath voltage, the secondary electrons are created on the surface of the electrode and more ionizations are observed in the sheath region. Meanwhile the electric field becomes stronger inside the sheath region. From Fig. 4(b), the secondary electron fluxes decrease quickly in the time range from 0 t/T to 0.45 t/T and increase quickly in the time range from 0.45 t/T to 1.0 t/T. The evolution of the secondary electron flux increases with increasing SEE coefficient. In the beginning and the end of one cycle, the evolution of the secondary electron flux increases largely with increasing SEE coefficient. From Fig. 4(c), the instantaneous ion flux increases quickly from 0 t/T to 0.45 t/T and decreases quickly from 0.45 t/T to 1.0 t/T. The evolution of the ion flux increases negatively with increasing SEE coefficient.
The electron pressure heating is mainly due to the electron diffusion.[23] Figure 5(a) shows that the cycle-averaged electron pressure heating increases with increasing SEE coefficient and mainly occurs in the sheath region. There are two negative peaks of the electron pressure heating in both sheath regions, and the peaks increase negatively as the SEE coefficient increases. According to Eq. (11), it is easy to find that electrons diffuse along the negative gradient direction of the electron density and move along the negative direction of the electric field. When these electrons reach the centre of discharge, where the electric field becomes low, they would not be cooled any more. In the bulk plasma region, the electron pressure heating is almost constant as the SEE coefficient increases.
The electron Ohmic heating is directly proportional to the square of the electric field, coupling with the electron density as can be seen in Eq. (12). The electric field is varied by adjusting the electrode gap, varying SEE coefficient or applying RF-voltage and frequency.[22] The electric field and electron density increase with temperature rising. It is suggested that the electron Ohmic heating becomes more effective with increasing SEE coefficient, which could theoretically result in escaping electron heating. From Fig. 5(b), it is worth noting that the higher electron Ohmic heating occurs in the sheath region and increases with increasing SEE coefficient. But the change of electron Ohmic heating is virtually constant in the plasma region. Also, the results show that the width of the sheath region shifts slightly toward the electrodes as the SEE coefficient increases. The electron Ohmic heating is large in the sheath region because a large number of electrons are accelerated by the electric field in the sheath region compared with in the bulk plasma region.
As described, most electron heating system configurations consists of two mechanisms: the electron pressure heating and the electron Ohmic heating as shown in Eq. (10). With increasing electron Ohmic heating, the chance of collisions between electrons and neutral particles increases, which leads to the growth of electron heating. Figure 55(c) shows the cycle-averaged electron heating. In the plasma region, the electron heating is positive and slightly increases as the SEE coefficient increases. In the two sheath regions, the electron heating is positive and increases as the SEE coefficient increases, except near the electrodes where the electron heating takes the negative values. The negative electron heating can be thought of as “diffusion cooling” since electrons diffuse against the local field in these regions.[23–25,31]
Figure 5(d) shows that the cycle-averaged ion heating changes with SEE coefficient in a range from 0.01 to 0.3. In the two sheath regions, the ion heating has a significant effect and increases as the SEE coefficient increases. But, in the bulk plasma region, the ion heating is almost zero. In the powered sheath region, the ion heating decreases quickly from the powered electrode. Also, in the grounded sheath region, the ion heating decreases quickly from the grounded electrode. The ion heating is positive where the location of ion-ion energy exchange is large because the energy exchange usually is faster than ion heating, and quite a quantity of charges can be produced there. There is a large interaction of the charge with plasma ion species and the interaction causes more ions to be heated. Also, ions are heated due to interactions between the electric field which collides with the sheath.[16–26]
The cycle-averaged electron energy losses are plotted in Fig. 6 as the SEE coefficient increases. The cycle-averaged electron energy loss occurs in the whole discharge region. The electron energy loss curves each have one peak with SEE coefficient of 0.01 and have two peaks when the SEE coefficient is bigger than 0.01 in the sheath region. It is known that the electrons lose energy as a result of the inelastic collsions between particles and electrons. The electrons gain energy from the electric field and lose energy by their colliding with neutral particles. The electron energy loss depends exponentially on electron temperature but linearly on electron density. It can be seen that as the SEE coefficient increases, the electron energy loss increases more in the bulk plasma region.
4. ConclusionsA 1D fluid model of capacitive RF argon discharge at low pressure is used to study the effect of the effective SEE coefficient γ on the glow discharge characteristics. The evolutions of the discharge parameters are investigated numerically. The results indicate that the electron density increases largely with electric field slightly increasing. Also the electron temperature increases in the two sheath regions and is almost constant in the bulk plasma region. The variations of electron current density between parallel-plate electrodes with the SEE coefficient increasing are presented and compared with each other. The evolutions of the secondary electron flux and the ion flux on the surface of the powered electrode are studied. With the SEE coefficient increasing, the negative maximum of the electron pressure heating in the sheath region increases while it is essentially unchanged in the center of the glow discharge. With SEE coefficient changing from 0.01 to 0.3, there appears a positive maximum of the electron Ohmic heating in each of the two sheath regions, but in the center of the glow discharge, it is almost constant. In addition, the electron heating increases, leading to a corresponding increase of the electron pressure heating and Ohmic heating in glow discharge. The results show that ion heating increases in the two sheath regions. The higher ionization rate and electron energy loss are found in the two sheath regions rather than in the bulk plasma region where it is found to be lower with all SEE coefficients. In general, these numerical results could provide valuable information about electrical glow discharge which otherwise could be determined only by experimental effort.